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£^ , Abstract 

The data from the QSO absorption lines indicating a nonzero time- variability of the 
. hne-structure constant has been re-analyzed on the basis of a "damped-oscillator" fit, 

£jy as motivated by the same type of behavior of a scalar field, dilaton, which mimics 

, a cosmological constant to understand the accelerating universe. We find nearly as 

1 good fit to the latest data as the simple weighted mean. In this way, we offer a way 

to fit the more stringent result from the Oklo phenomenon, as well. 

CM ■ 1 Introduction 

> 

In this article we suggest a possible way to understand the nonzero time- variability of the fine-structure 

stringent constraint from the Oklo phenomenon [3-5]. Interesting enough this consistency is derived 
from the suspected behavior of the scalar field required to fit the accelerating universe [6], another 
striking finding in today's cosmology. The rest of this Section will be devoted to a brief account on 
the background as well as the summary of the result. 



Of 



^ ■ Attempts to probe Aa for the difference from today's value of the fine-structure constant a based 

on the calculated relativistic corrections of many transitions is now summarized [2] to give Aa/a = 
(-0.54 ± 0.12) x 10~ 5 for the redshift range 0.2 - 3.7, or the fractional look-back time 0.2 - 0.9. 
The estimate according to the weighted mean of the 128 data points, which we accept in the present 
analysis, may be referred to as the 1-parameter (horizontal straight-line) fit, for the convenience of 
comparison with our own fit in terms of a "damped-oscillator" model with three parameters, as will 
be explained shortly. 

The 1-parameter fit in this sense appears to be confronted by the constraint smaller by 2 or 3 orders 
of magnitude as was established by the measured isotopic abundance of 149 Sm in the remnant of Oklo 
natural reactors. This result, which, despite apparent complications of nuclear physics, is relatively free 
from uncertainties thanks to an exceptionally low-lying resonance, is confirmed by careful re-analyses 
of Ref. [3], yielding either of the upper bounds |Aa/a| <,10 -7 [4] or <10~ 8 [5], or the nonzero value 
~ 10~ 7 [5] still not completely ruled out. The small ratio of the resonance energy 97.3 meV to ~MeV 
that characterizes typical energy scale of nuclear phenomena makes it also plausible to derive the result 
on Aa/a without entering into details of the strong-interaction coupling constant [5]. 

No immediate conflict is encountered, however, because the fractional look-back time around 0.14 
for the Oklo phenomenon is outside the range for the QSO observation. The small time-separation 
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between them still seems to indicate a rather sharp time-dependence of a, obviously requiring a the- 
oretical interpretation. Note also that it imposes a much more stringent constraint than the results 
from such low-z phenomena, like Hi 21cm and molecular rotational transitions [7]. 

It was argued that the issue can be related to the accelerating universe which is accounted for by 
a cosmological constant with ~ 0.7 [6] by allowing the time-dependent speed of light [8]. Also 
shown was the potential importance of spatial inhomogcncity of time variation of a, concluding that 
the phenomena taking place on the Earth's surface might be affected by the local matter density more 
than those in the galactic clouds [9] . We pursue here instead a more conventional approach by making 
use of somewhat nontrivial time-dependence suggested by the behavior of the scalar field supposed 
to result in an effective cosmological constant. Possible sharp change of the scalar field has also been 
discussed in connection either with the running nature of its coupling constant [10], or with a classical 
field, like ours but motivated with different potentials [11,12]. 

Our detailed analysis of the accelerating universe in terms of the scalar-tensor theory [13,14] in- 
dicates a particular type of behavior of the scalar field, sometimes called a dilaton a, which stays 
nearly constant that mimics a cosmological constant. A closer look shows further that this "plateau" 
is superimposed with a small vibration, reminiscent of a damped oscillation as a function of the cosmo- 
logical time, as will be explained briefly in Section 2. We combined this with a relation Aa/a oc Act, 
which we derive based on a specific model, but might be viewed as an extension of the previous results 
[15,16]. We try to see if this oscillation pattern can be fitted to the observed Aa/a with a zero of the 
oscillation identified with the small value constrained by the Oklo phenomenon [17]. 

In view of still unavoidable uncertainties in selecting a unique solution of the cosmological equations, 
we attempt a phenomenological approach by assuming the damped oscillation behavior of Aa/a with 
three arbitrary coefficients which arc determined to fit the QSO observation. We show that this 3- 
parameter fit to the QSO result can be nearly as good as the 1-parameter fit in [1,2], fitting the Oklo 
data as well. Notice that the 1-parameter fit, if extended to the Oklo time, would yield an unacceptably 
large value of x 2 ■ 

One might still wonder why we are content with X 2 e d = 1-09 for the former with more adjustable 
parameters is not significantly smaller than x 2 e d = 1-06 for the latter. In reply, we are content, at this 
moment, with the result that our fit motivated by a theoretical prejudice survives a realistic test. This 
"success" will be followed by future attempts to scrutinize the solutions for the accelerating universe. 
Admitting that the result is far from final, we still believe it to serve as a useful first step toward better 
understanding the issue in a wider perspective. 

In Section 2, we begin with describing the current theoretical model following [13,14]. Section 3 
offers our 3-parameter fit and presenting the result, to be compared with the starting cosmological 
solution. The final Section 4 is devoted to concluding remarks, including a discussion on the constraint 
proposed recently on the Re-Os decay [18]. 
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2 Theory of the cosmological constant and the time-dependent 
fine-structure constant 



Today's cosmological constant problem has two faces [14]: Why is A bs so small compared with the 
theoretical expectation by 120 orders of magnitude? Why is it still nonzero? 

The scalar-tensor theory provides an acceptable answer to the first question. The theory seems to 
deserve two renewed interests; possible origin of the scalar field in the dilaton of string theory, and 
a way to implement the scenario of a decaying cosmological constant. One of the crucial theoretical 
ingredients is in choosing a correct conformal frame (CF). String theory is formulated naturally in the 
J(ordan) CF, featuring a "nonminimal coupling" of the dilaton field (though in higher dimensions). 
We also accept that, almost in any version of unification models, a cosmological "constant" A is present 
in the J frame, namely in the "string CF." 

A "physical" CF is to be chosen, on the other hand, according to the two requirements; i) keeping 
the solution of the cosmological equations with A included from deviating too much off the standard 
cosmology, ii) resulting in the masses of matter fields time-independent. The second condition is 
imposed because our clocks and meter sticks are made basically of the particle masses. A solution to 
this non-trivial problem is offered by assuming a "scale-invariance" in 4 dimensions except for the A 
term, finding the E(instein) CF with the Einstein-Hubert term as physical. Particle masses emerge 
nonzero as the invariance is broken spontaneously. The A term in the J frame is now converted to an 
exponential potential of the scalar field a. Unlike in the more phenomenological approach with the 
quintessence models, no inverse-power potential is expected. The energy density of a, dark energy, 
falls off like t~ 2 , thus implementing the scenario of a decaying (effective) cosmological constant [19]. 

The scalar field a remains decoupled from matter in the classical limit, but is brought to matter- 
coupling due to the interactions among matter fields, a la quantum anomaly. According to the simpli- 
fied QCD calculation, the coupling constant to the quark is obtained; g aq — ((5a s /n) rts 0.3£, where 
( is a constant of the order one, a s w 0.2 being the QCD analog of the fine-structure constant, while 
Mp 1 = 8ttG has been suppressed. The fact that the matter coupling depends on the strong- interaction 
constant signals composition-dependence of the force mediated by a, leading to WEP violation. 

The coupling constant to nucleons is also computed; g a ^ ~ 0.019£. The relatively small value is 
a consequence of the smallness of the quark-mass component in the nucleon mass. The a eventually 
acquires a nonzero mass, as a pseudo Nambu-Goldstone boson, resulting in a medium-ranged force. An 
exact estimate of the force-range A is not easy, but is certainly shorter than distances that characterize 
solar-system experiments, which therefore no longer constrain the parameters of the a force. 

The parameter a 5 in the non-Newtonian potential between two nucleons —{Gm^/r){\ + a<se~ r / x ) 
[20] is given by a$ — 2g^ N w 0.72 x 10~ 3 £ 2 , being consistent with the composition-independent 
experiments roughly for A<lm [21] including the recent submillimeter experiment [22]. Results of 
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composition-dependent experiments may be analyzed by replacing 05 — ► a^q^qj with the "charge" g» for 
the i-th substance. For a free-fall experiment on the Earth, we find the relative acceleration difference 
(Ag)ij t E/g ~ Q!5(Ag)jjA x 4 x 10 -8 , where A is measured in units of meters. In view of still present 
uncertainties in (Aq)ij, we borrow tentatively from the vector fifth-force model, \(Aq)ij\ ~2x 10~ 3 
[21]. The available observational upper bounds \(Ag)ij^/g\ ~ 1CT 12 is respected if A<10m, thus giving 
an idea how large the expected WEP violation can be. 

The second face, or question, requires some deviation from the smooth behavior of the energy 
density p a . As one of the simplest approaches along this direction, we introduced another scalar field, 
X, assumed to have a specific interaction with a. A typical cosmological solution for k = is shown in 
Fig. 1 around the present epoch (taken from Fig. 5.8 of [14]). The energy density p s of the dark energy 
<t, x behaves like a plateau, acting as a cosmological constant. As a consequence of the dynamics of the 
scalar fields, the plateau comes to a crossing with the ordinary matter density. Nearly in coincidence 
with this, the scale factor exhibits a "mini-inflation," that fits the observed "acceleration," also with 
concomitant sudden changes of a and \- 




Figure 1: An example of the Friedmann solution 
with k — of the cosmological equations of the 
two-scalar model, taken from Fig. 5.8 of [14]. On 
the horizontal axis, logi ss 60 corresponds to the 
present epoch in the (reduced) Planckian unit sys- 
tem of c = h = M P (= (8ttG)~ 1/2 ) = 1. In the lower 
panel, p and p s for the ordinary matter density and 
the energy density of the dark energy, namely for a 
and x 5 respectively are shown. 
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Figure 2: Magnified view of a (solid) and 0.02x + 
44.25 (dashed) vs \ogt near the present epoch. Note 
that before magnification, this is part of the plot 
(previous figure) in which a, after the period of 
plateau behavior resulting in a "mini-inflation" of 
the scale factor, shoots up "suddenly." 



Particularly interesting is found if the apparently simple jump of a, as shown in Fig. 1, is vertically 
magnified greatly, by 330 times, for example, as shown in Fig. 2. This unexpected damped-oscillator- 
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like behavior plays a key role in bringing about an accelerating universe. The oscillation is "small," 
however. We recognize it only after the plot is magnified considerably than what we need to analyze how 
the cosmological acceleration takes place. Different cosmological solutions leading to nearly identical 
results may be differentiated only by looking at other phenomena, including Aa/a, among others. 
We derived the relation (6.194) of [14] through the QED calculation; 



result only as an example of the expected relation on the idea that a possible variation of the coupling 
constant is due to the changing scalar field, as we may obtain directly from string theory, though 
uncertainties seem unavoidable for deriving results applicable directly to 4 dimensions. In the following 
we take a more phenomenological attitude that the pattern of a(t) is described by a simplified damped 
oscillation with the parameters determined to fit the observation of Aa/a. 

3 Analysis of the data 

We begin with Fig. 2 showing a and % as functions of logt, which is a natural choice of the time 
variable in solving the cosmological equation. The function o is readily converted to that of the 
fractional look-back time u = (to — t)/to, with to the present age of the universe [1,2], essentially in 
the backward time direction. From a phenomenological view suggested above, we assume the purely 
damped oscillation behavior for y = (Aa/a) x 10 5 ; 



where u ki is chosen to be 0.142 corresponding to the mean (1.95 ± 0.05)Gy of the time of the Oklo 
phenomenon. The natural reactors are believed to have lasted somewhat shorter than 10 6 years, 
corresponding to a width in the u-direction as small as ,<10~ 4 . According to a positive (negative) a, 
y(u) increases (decreases) at u = u kio- We should choose b > 0. 

We notice that choosing (2) results in y(0) = asm (—2t:u ^\ /T) which may not vanish in general, 
in contradiction with the natural condition y(0) = to be obeyed by definition. As the simplest 
remedy we may add an offset constant to the right-hand side of (2), but calling for another parameter. 
Improving the analysis in this direction will be made in the future, however, by emphasizing the 
behavior in the ranges close immediately to the QSO range, for the moment. 

The data set to be fitted is taken from the "fiducial sample" (Fig. 8) of [2] of 128 points. A 
1-parameter fit results in 



with Xrcd = x/127 = 1.06, which would have been 1.24 if we were to try to fit the data by the null 
result y = 0. We emphasize that the error bar at the Oklo point is invisibly small (<,10 -2 ) compared 




Aa ^ a A A 

= Z— (Aa : 

a 2ir 

larged fields tha 




(2) 



y = -0.54 ±0.12 



(3) 
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with other error bars. It then follows that the 1-parameter fit if extended to u kio would have resulted 
in x^ ed as large as the order of >,10 4 . In order for the result to be sensible, we must bend the straight 
line (3) to pass the point Uokio on the u-axis, but hopefully with a theoretical rationale. 



Solution 


b 


a 


T 


Xrod 


A+ 


2.4 


0.151 


0.714 


1.09 


B+ 


0.0* 


0.350 


0.275 


1.22 


c+ 


3.2 


0.041 


0.182 


1.22 


A_ 


4.0* 


-0.035 


0.451 


1.20 


B_ 


2.2 


-0.085 


0.218 


1.21 


C_ 


3.0 


-0.037 


0.158 


1.23 



Table 1: Six solutions for minimized \ 2 y an< i their 
parameters. The asterisks mean the end-point min- 
ima instead of the genuine local minima. 



b 


a 


T 


x 2 


^rcd 


0.0 


0.728 


0.716 


136.83 


1.10 


2.4 


0.151 


0.714 


136.19 


1.09 


4.0 


0.051 


0.710 


136.46 


1.09 



Table 2: Parameters for the local minima for the 
solution A+, corresponding to the three curves in 
Figs. 3 and 4. In the second line for b — 2.4, the val- 
ues are for the 3-dimensional minimum of x 2 5 while 
other two lines for b — 0.0 and b — 4.0 represent 
only the 2-dimensional minima listed for comparing 
X 2 with the true minimum for b — 2.4. 



Using the behavior of the scalar field for the accelerating universe as a guide, we conveniently 
consider the range < b < 4, inside of which we search for sets of parameters a, b, T that minimize \ 2 - 
The Oklo point is not included in x 2 because choosing it o ki is meant to fit the data already. We have 
found six solutions, with the parameters and x 2 c d = X 2 /125 summarized in Table 1. Remarkably, only 
the solution A + offers x 2 c d = 1-09 which is small enough nearly comparable with the value 1.06 for the 

1- parameter fit (3), while others share larger values ~ 1.2 close to 1.24 for the null result mentioned 
above. There are yet other solutions with even larger \ 2 ■ In what follows we focus on A + . 

We consider x 2 in 3-dimensional space of a,b,T. For each solution there is a 3-dimensional local 
minimum, around which there occurs a region of confidence level of 68%, for example. For A + , the 
true minimum occurs for b = 2.4, a = 0.151, T — 0.714, as given by the first line of Table 1, which 
are also reproduced in the second line of Table 2. The surface which surrounds the confidence volume 
may be represented by the three cross sections each of which is a contour in the 2-dimensional a — T 
plane with constant b, as shown by Fig. 3 and Table 2. The complete 3-dimensional minimum lies 
on the plane of b = 2.4, as marked by a cross inside a contour labeled by 2.4 in Fig. 3, while the 

2- dimensional minima are also shown as dots inside the other two contours for b = 0.0 and b = 4.0. 
The corresponding values of x 2 are given in Table 2. The largest contour occurs for b = 0.0, though 
the corresponding minimum x 2 is larger than that of b = 2.4. In this connection we further point out 
that b = 0.0 is likely outside the practical interest, because we expect that the present epoch is well 
before the time of sudden changes of the scalar fields, as shown in Fig. 2. 

Figure 4 shows how the (solid) curve of A + fits the data. Also shown are the curves of neighboring 
fits corresponding to b = 0.0 (dotted) and b = 4.0 (dashed), which nearly overlap that of b = 2.4 in the 
range of the QSO observation. There is an indication, however, that the larger b makes it the easier 
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to fit the Oklo result. In fact choosing b — 8.0(a = 0.00297, T — 0.712) gives the curve to be within 
±0.001, hence satisfying the "stronger" Oklo constraint |Aa/a|;< 10~ 8 in the range u = 0.142 ± 0.02. 
This is made, however, with the cost of somewhat larger xj: ed = 1.11, because the QSO data exhibits 
quite "flat" distribution, instead of a sharp change implied by a large b. In fact this flatness makes it 
unlikely that the smaller value at the Oklo point is largely due to the damping effect, e bu , even taking 
aside that a large b is not favored by the cosmological behavior. 




0.2 0.4 0.6 0.8 1 




0.4 0.6 0.E 

Fractional look-back time u 



Figure 3: Cross sections by a ~ T planes, labeled by 
the values of b, of the 68% confidence level of the 
solution A+ . The true 3-dimensional minimum of \ 2 
occurs for the point denoted by a cross inside the 
contour labeled by 2.4, as also given in Table 2. The 
points denoted by dots inside other two contours rep- 
resent only the 2-dimensional minima for each b. 



Figure 4: Fit of the solution A+ (solid). Similar 
curves (dotted for b — 0.0, dashed for b = 4.0) are 
also shown corresponding to Table 2. In the upper 
panel, the curves for u < 0.2 have been 10 times 
magnified. In the lower panel (binned data), the 
(dot-dashed) horizontal straight line is added cor- 
responding to the 1-parameter fit. 



4 Concluding remarks 

We have shown that the QSO data can be fitted by a damped-oscillator fit nearly equally well as 
the 1-parameter fit. Also we can do this by fitting the Oklo result as well. We have, however, some 
unsolved problems. 

The behavior as shown in Fig. 1 indicates b « 2.5 and T « 0.22, compared with b = 2.4 and 
T = 0.714 for the fit A + ; a near agreement for b but the latter T seems somewhat larger than the 
former. Smaller T tends to give larger Xrcd as we infer from Table 1. Also the maximum of |A<r| 
around logi w 59.6 for u w 0.7 is ~ 0.1, which should be compared with |t/| max ~ 1 indicated in Fig. 4. 
In the cosmological solution in [14] we used the value Q — 1.58. Combining these with Z = 5 from the 
standard gauge theory, we find that Aa/a in A + is smaller than the right-hand side of (1) by about 
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2 orders of magnitude. 

We add a comment on the /3-decay of 187 Re, which is recently claimed to require \Aa/ce\<3 x 10~ 7 
for the last 4.5 Gys, or u < mr c w 0.33 [18]. The solution A + and its neighbors shown in Table 2 give 
the values that fail to comply this constraint. Even the choice of b — 8.0 as before gives 4.3 x 10~ 7 , 
which comes only close to the required bound. We may reverse the analysis by first identifying the 
rhenium constraint with a zero of y, and finding (Aa/a)oido ~ 1-0 x 10~ 7 . This might be accepted as 
the largest possible deviation of the Oklo constraint, with x^cd = 1-09, but the required b — 8.0 (also 
with a = —0.00366, T = 1.052) is again quite outside the range allowed from cosmology. 

It appears that the rhenium constraint can be made only barely consistent with the flat distribution 
of the QSO result, at the best. The same can be a problem for other non-oscillatory behaviors as well. 
The function used in [10], for example, features a sharp decrease of |Aa/a| from z = 2 to z = 0.45 
for the rhenium time to account for the small value. This entails, in turn, a sharp increase toward the 
high-z end of the QSO range, agreeing only marginally with the data including not only the average 
value but also more details, which can be exploited by a x 2 analysis, for example. A conclusion of the 
similar nature will follow also from Fig. 3 of [11] and Fig. 13 of [12], for which the estimated x^ed 
appear to remain within an acceptable range, but |Aa/a| calculated at the rhenium time turn out 
somewhat larger than the target value. 

In this connection we point out that estimating the parameters of Re-Os isochron accurately is 
primarily for obtaining the time-averaged decay-rate [23] eventually to determine the age of meteorites. 
Compared with the Oklo phenomenon and the QSO absorption lines, the possible time- variation of the 
decay-rate, namely of a, at certain time has been inferred only by the argument which is less direct 
and less compelling, as will be elaborated separately [24]. 

As a unique feature of the current model, we point out that the behavior of a toward the present 
epoch depends sensitively on the status at much earlier times. Figure 5.14 of [14] provides examples 
to show how small and smooth changes of the initial value may affect whether an oscillation occurs 
or not, and if it does, how much the amplitude and the frequency vary. The example shown in Fig. 
2 turns out to have a relatively large amplitude. A much smaller amplitude can be chosen without 
noticeably affecting the way the universe is accelerated. Also the interaction between a and x was 
assumed to be the simplest ((5.58) of [14]) to implement the mechanism for the acceleration. One 
might even modify this interaction in such a way that fitting the Oklo data as well as meeting the 
condition y(0) = can be made more natural. These details will be the subject of future studies. We 
re-emphasize that measuring Ace /a, even if it might finally turn out to be smaller, will constrain the 
behavior of the scalar field rather independently of its role for the cosmological acceleration. 

It seems nevertheless worth emphasizing that the size of the prediction (1) is not very much away 
from the observation, allowing us to start with this relation as a reasonable basis. The issue is more 
acute when we notice that the QSO result implies a/ a ~ 10 _16 y _1 , and its magnitude, several orders 
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smaller than t$ ~ 10 _10 y _1 , has never been fully understood from a theoretical point of view. 

I thank Michael Murphy for providing the QSO data before the publication and reading the 
manuscript, and also Takashi Ishikawa for the helpful discussion on the error estimates. My thanks 
are also due to Akira Iwamoto for his crucial comments on the rhenium decay. 
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